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CLAIMS 

What is claimed is: 

1. (Cancelled) 

2. A method for computing the conformal structure of a surface, 
the method comprising: 

receiving a first mesh representation M of said surface; 

conformally mapping said mesh representation to a canonical 
parameter domain forming a first mapped surfaces; and 

computing a conformal parameterization of said first mapped 
surface. 

3. The method of claim 1 wherein the step of computing said 
conformal parameterization includes: 

computing a homology basis {yi, yzt - / Y2g} of said mesh 
representation, wherein computing said homology basis includes 
the steps of: 

computing a set of boundary matrices di, dz for said 
mesh representation^ respectively; 

forming a matrix D, wherein D = dzdi + didi; 

computing the Smith normal for said matrix D; and 

computing the eigenvectors of D corresponding to zero 
eigenvalues that are equal to {yi, ya, ... , yig) , and wherein {yi, y2, 
... , y2g} are a homology basis of said mesh representation. 

4 . The method of claim 2 wherein the step of computing the 
conformal parameterization includes the step of computing a set 
of harmonic l-form basis, wherein the step of computing said set 
of harmonic 1-form basis includes the steps of: 
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calculating the value of ci^ = -yi • Yjf i/j = 1, 
2, 2g wherein ji, Yj is the homology basis of the mesh 

representation; and 

solving the linear system bOi = 0, Atoi = 0, <G)i, yj> 

5 = -Yi • Yj for coi, wherein {qi, C02, ... , 0)2g} is the desired harmonic 
1-form basis of said mesh representation. 

5. The method of claim 3 wherein the step of computing the 
conformal parameterization further includes the steps of: 
10 computing the fundamental domain Dm of the mesh 

representation, wherein the step of computing the fundamental 
domain Dm of. the mesh representation includes the steps of: 
selecting an arbitrary face, fo e M; 
setting Dm to fo; 
15 setting SDm = dfo; 

placing all neighboring faces of fo that share an edge with 
fo in a queue, Q; 

while Q is not empty: 

reading the first face in Q; 
20 setting af = eo + ei + e2; 

setting Dm = Dm U f; 

finding the first ei e df such that -ei e aOn; 
replacing -ei in ^Dm with (ed+D, e(i+2)}; 

putting in Q all the neighboring faces that share an 
25 edge with f in the mesh representation and that are not in Dm or 
Q; and 

removing all adjacent oriented edges in SDm that are opposite 
each other in sign. 
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6. The method of claim 4 wherein the step of calculating said 
conformal parameterization of said first mapped surface further 
includes the steps of: 

recording a first path 5 from root vertex vo to u by 

5 traversing all the vertices u e Dm; 

computing the integration of (|)(u) = <Ci>, yu>; and 

providing as an output ^(u) as the conformal coordinates of 

u. 

10 7. The method of claim 5, wherein the step of calculating said 
conformal parameterization of said first mapped surface further 
includes the step of calculating the cohomology basis, wherein 
the step of calculating the cohomology basis includes the steps 
of: 

15 setting C0i(ei)=l and a)i(ej)=0 for any edge e € T; 

ordering Dm such that Dm = {fi, ±2, r fn); 
reversing the order of Dm to {fn, fn-i,f fil; 
while Dm is not empty: 

retrieving the first face f of Dm; 
20 removing f from Dm, 9f = eo + ei + 62; 

divide {ek} into two sets, r = {eGSf l-eeaD^}, 

choosing the value of (Oi(e,^),e,j €11 arbitrarily, such that 
^^^^a)i(e)=-]^^^^Oi(e), and if 11 is empty, then the right hand side 

25 is equal to zero; 

updating the boundary of Dm/ let 9Dj^ =9Dj^ +5f ; and 
computing the dual of the homology basis, {yi,Y2»— >Y2g}' 
the linear transform {to,,a)2,...,co2g} such that <yi,G)j >=-yj yj . 
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8. The method of claim 1, wherein the step of conformally 
mapping and calculating said conformal parameterization of said 
mesh representation includes the step of calculating the 
holomorophic 1-form, wherein the step of calculating the 
5 holomorphic 1-form includes the steps of: 



9. The method of claim 7 wherein the step of doubling M to form 

M includes the steps of: 



reversing the orientation of -M; 

finding for each boundary vertex u e 3M, a unique 
corresponding boundary vertex -u € 9 -M; 

finding for any edge on e e 9M a unique boundary edge -e e d 
20 -M; and 

gluing M and -M such that the corresponding vertices and 
edges are identical, wherein the resulting mesh is the doubling 



25 10. The method of claim 7 wherein the step of calculating said 
conformal parameterization of said mesh representation, M, 
further includes the step of calculating a conformal structure, 
wherein the step of calculating the conformal structure includes 
the steps of: 

30 a) providing tlie holomorphic 1-form basis {r/+V-Tr,*|; 



10 




15 



copying M, denoted as -M; 



M. 
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b) computing a partition (Ui), such that M c Ui, Ui being 
simply connected; 

c) for each Ui, selecting a holomorphic 1-form base 

(Dj + V-Io)] ; 

5 d) integrating the holomorphic 1-form on Ui, denoting the. 

mapping as zi, and in the event that there are zero points, 
subdivide Ui and repeat the integration; and 

e) outputting the conformal structure { (Ui, zi) } . 

10 11. The method of claim 4 wherein the step of calculating said 
conformal parameterization of said mesh representation further 
includes the step of calculating the period matrices, wherein the 
step of calculating the period matrices includes the steps of: 

calculating the elements of a matrix C as cij = <yi, coj>; 

15 calculating the elements of a matrix S as sij = <Yi, 

C0j*>; and 

calculating the matrix R as the solution of CR = S where R 

satisfies = -I, where I is the identity matrix. 

12. A method for classifying a surface, the method comprising 
the steps of: 

receiving a mesh representation of the surface; 
computing a period matrix R corresponding to said mesh 
representation of the surface; and 
storing said period matrix R. 

13. The method of claim 11 wherein the step of computing said 
period matrix R includes the steps of: 

computing a homology basis {yi, 72, ... / 729} of said mesh 
30 representation, wherein computing said homology basis includes 
the steps of: 



20 



25 
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computing a boundary matrix di, d2 for said mesh 
representations- 
forming a matrix wherein D = 3232^ + di'^di; 
computing the Smith normal for said matrix D; and 
5 computing the eigenvectors of D corresponding to zero 

eigenvalues that are equal to {yi, 72/ ... / Jig) wherein {yi, 72, ... , 
y2g} are a homology basis of said mesh representation. 

14. The method of claim 12 wherein the step of calculatihg said 
10 period matrix R further includes the step of computing a set of 
harmonic 1-form basis, wherein the step of computing said set of 
harmonic 1-form basis includes the steps of: 

calculating the value of ci^ = -yi • yj, i, j = 1, 2, 

2g, wherein ji, yj are the homology basis of the mesh 
15 representation; and 

solving the linear system 5(0i = 0, Acoi = 0, <(^i, yj> 

= -yi • yj for coi, wherein {coi, 0)2/ - / <^2q) is the desired harmonic 
1-form basis of the mesh representation. 

20 .15. The method of claim 13 wherein the step of computing said 
period matrix R further includes: 

calculating the elements of a matrix C as cij = <yi, C0j>; 

calculating the elements of a matrix S as Sij = <yi, 
(i)j*>; and 

25 calculating the period matrix R as the solution of CR = S 

where R satisfies R^ = -I, where I is the identity matrix. 

16. A method for matching first and second surfaces, the method 
comprising: 

30 receiving first and second mesh representations of said 

first and second surfaces, respectively; 
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conformally mapping both of said first and second mesh 
representations to a first canonical parameter domain, thereby 
forming first and second mapped surfaces; 

computing first and second conformal parameterizations of 
5 said first and second mapped surfaces, respectively; 

computing first and second level sets of Gaussian curvature 
and mean curvature for said first and second mapped surfaces, 
respectively, wherein said first and second level sets of 
Gaussian curvature and mean curvature are a function of said 
10 first and second conformal parameterizations, respectively; and 

comparing said first and second level sets of Gaussian 
curvature and mean curvature and in the event that the comparison 
exceeds a predetermined threshold, declare a match between said 
first and second surfaces, otherwise declare a mismatch between 
15 said first and second surfaces. 

17. The method of claim 15 wherein the step of computing said 
first and second conformal parameterizations includes the step 
of: 

20 computing a homology basis {yn, Yi2/ ... / Yi2g} and {yai, 722^ 

/ Y22g} of said first and second mesh representations, 
respectively, wherein computing said homology basis includes the 
steps of: 

computing first and second sets of first and second 
25 boundary matrices Sn, diz and 821, 822 for said first and second 
mesh representations, respectively; 

forming first and second matrices Da and D2, 
respectively, wherein Di = 5i2<9i2^ dii^dn and D2 = ^22^22 + 

<52i^^2i; 

30 computing the Smith normal for said matrices Di and D2; 

and 
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computing the eigenvectors of Dl and D2 corresponding 
to zero eigenvalues that are equal to - {yn, yi2, ... , yi2g} and 

{y2i/ y22/ ... . 722g) wherein {yxi, yx2/ ... , yx2g} are a 

homology basis of said first and second mesh representations, 
5 respectively. 

18. The method of claim 16 wherein the step of computing said 
first and second conformal parameterizations further includes the 
step of computing a first and second set. of harmonic 1-form 

10 basis, wherein the step of computing said first and second set of 
harmonic 1-form basis includes the steps of: 

calculating the value of Cn^ = -yu • yij, ij = 1/2, ... 

, 2g and Czi^ = -yzi • y2j, i, j = 1. 2, ... , 2g wherein yn, yij and y2i, 
y2j are the homology basis of the first and second mesh 
15 representations, respectively; 

solving the linear system Scoii = 0, AcDii = 0, <coii, yij> = -yn • 

yij for Oii, wherein {con, 0)12, ... , 00129} desired harmonic In- 

form basis of the first mesh representation; and 

solving the linear system 5co2i = 0, Aco2i=0, <co2i, y2j> = -y2i • 
20 y2j for Q2if wherein {CO21, 0)22/ r ^22q} is the desired harmonic 1- 
form basis of the second mesh representation. 

19, The method of claim 17 wherein the step of computing said 
first and second conformal parameterizations further includes the 

25 step of computing the fundamental domain Dim and D2M of the first 
and second mesh representations respectively, wherein the step of 
computing the fundamental domain .Dim and D2M of the first and 
second mesh representations respectively includes the steps of: 

selecting an arbitrary face, fio g Mi and f2o e M2, and set 
30 Dim and D2M to fio and f2o respectively; 

setting ^Dim = 5fio and aD2m = 5f2o; 
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placing all neighboring faces of fio and ±20 that share an 
edge with fio and fao respectively in a queue, Qi and Q2, 
respectively; 

while Qi and Q2 are not empty: 

reading the first face in Qi and Q2; 
setting dfi = eio + en + 612; 
setting Dim = Dim ^ fi; 

finding the first en € 3fi such that -en e aOiM; 

replacing -en and --e2i in dDm and aD2M/ respectively, 
with {ei(i+i), ei{i+2)} and {e2(i+i)/ 62(1+2))/ respectively; 

putting all the neighboring faces that share an edge 
with fi and f2 in the first and second mesh representations, 
respectively, and that are not in Dim and Qi and D2M and Q2, 
respectively, into Qi and Q2, respectively; and 

removing all adjacent oriented edges in dDm and (9D2m that are 
opposite each other in sign. 

20. The method of claim 15 wherein said first and second 
surfaces are human faces. 

21. The method of claim 15 wherein said canonical parameter 
domain is a disk. 

22. A method for recognizing a surface represented by a mesh 
representation, M, the miithod comprising the steps of: 

receiving the mesh representation of the surface; 
removing at least one feature point from the mesh 
representation; 

doubling the mesh representation, M; . 

computing the period matrix, R, for the mesh representation, 
M, having at least one feature point removed; 
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comparing the period matrix, R, corresponding to the mesh 
representation . with at least one standard period matrix 
corresponding to a standard surface; 

in the event that the period matrix, R, corresponding to the 

5 mesh representation and the at least one standard period matrix 
corresponding to a standard surface are within a predetermined 
factor of one another, providing indicia that the surface 
represented by the mesh representation is recognized as being 
substantially similar to the at least one standard surface. 

10 

23. The method of claim 21, wherein the step of doubling the 
mesh representation, M, includes the steps of: 

copying M, denoted as -Mr- 
reversing the orientation of -M; 
15 finding for each boundary vertex u e 5M a unique 

corresponding boundary vertex -u e d -M; 

finding for any edge on e e 9M a unique boundary edge 

-e e 9 -M; and 

gluing M and -M such that the corresponding vertices and 
20 edges are identical, wherein the resulting mesh is the doubling 

M . 

24. The method of claim 21, wherein the step of computing the 
period matrix R includes the steps of: 

25 computing a homology basis {yi, y2, ... r Y2g} of said mesh 

representations, wherein computing said homology basis includes 
the steps of: 

computing a set of boundary matrices Si, dz for said 
mesh representations; 
30 forming a matrix D, wherein D = 3232^ + di^dif 

computing* the Smith normal for said matrix D; and 
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computing the eigenvectors of D corresponding to zero 
eigenvalues that are equal to {yi, 72/ f y2g} wherein {yi, y2, ... , 
y2g} are a homology basis of said mesh representation. 



5 25. The method of claim 23 wherein the step of calculating said 
period matrix R further includes the step of computing a set of 
harmonic 1-form basis, wherein the step of computing said set of 
harmonic 1-form basis includes the steps of: 

calculating the value of Ci^ = -yi • yj, i, j = ' 1, 2, 

10 ... , 2g, wherein yi, y^ are the homology basis of the mesh 
representation ; and 

solving the linear system Scoi = 0, Aoi = 0, <a)i, yj> = -yi • yj 
for coi, wherein {coi, (02, ... , 0)29} is the desired harmonic 1-form 
basis of the mesh representation. 

15 

26. The method of claim 24 wherein the step of computing said 
period matrix R further includes: 

calculating the elements of a matrix C as cij = <yi, coj>; 

calculating the elements of a matrix S as Sij = <yi, 
20 (Dj*>; and 

calculating the matrix R as the solution of CR = S where R 
satisfies R^ = -I, where I is the identity matrix, 

27. A method for recognizing a surface represented by a mesh 
25 representation, M, the method comprising the steps of: 

receiving the mesh representation of the surface; 
removing substantially all feature points from the mesh 
representation; 

selecting an arbitrary point on the surface of the mesh 
30 representation, M; 
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selecting an arbitrary orbit along the surface of the mesh 
representation, M, for the arbitrary point to follow; 

moving the arbitrary point along the arbitrary orbits- 
removing the arbitrary point from the mesh representation, 
5 M, at at least one point along the arbitrary orbit; 

doubling the mesh representation, M, having all feature 
points removed and the arbitrary point removed; 

computing at least one period matrix, R, for the mesh 
representation having the arbitrary point removed; 
10 comparing the at least one period matrix, R, corresponding 

to the mesh representation with at least one standard period 
matrix corresponding to a standard surface; and 

in the event that the period matrix, R, corresponding to the 
mesh representation and the at least one standard period matrix 
15 corresponding to a standard surface are within a predetermined 
factor of one another, providing indicia that the surface 
represented by the mesh representation is recognized as being 
substantially similar to the at least one standard surface, 

20 28. The method of claim 26, wherein the step of doubling the 
mesh representation, M, includes the steps of: 
copying M, denoted as -M; 
reversing the orientation of -M; 

finding for each boundary vertex u e 5M, a unique 
25 corresponding boundary vertex -u e 5 -M; 

finding for any edge on e € 9M a unique boundary edge 
-e e 5 -M; and 

gluing M and -M such that the corresponding vertices and 
edges are identical, wherein the resulting mesh is the doubling 

30 M . 
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29. The method of claim 26, wherein the step of computing the 
period matrix R includes the steps of: 

computing a homology basis {yi, 72/ ... , Jiq) of said mesh 
representation, wherein computing said homology basis includes 
5 the steps of: 

computing a boundary matrix di, 82 for said mesh 
representation; 

forming a matrix D, wherein D = 3232^ + di^di^ 
computing the Smith normal for said matrix D; and 
10 computing the eigenvectors of D corresponding to zero 

eigenvalues that are equal to {yi, y2, ... , y2g} wherein {yi, ya, ... , 
y2g} are a homology basis of said mesh representation. 

30. The method of claim 26 wherein the step of calculating said 
15 period matrix R further includes the step of computing a set of 

harmonic 1-form basis, wherein the step of computing said set of 
harmonic 1-form basis includes the steps of: 



... , 2g, wherein yi, yj are the homology basis of the mesh 
20 representation; and 

solving the linear system 6c0i = 0, AcOi = 0, <G)i, yj> = -yi • yj 
for coi, wherein {qi, 0)2, ... , co2g} is the desired harmonic 1-form 
basis of the mesh representation. 

25 31. The method of claim 29 wherein the step of computing said 
period matrix R further includes: 



calculating the value of ci^ = -yi • yj, if j = 



Ir 2, 



calculating the elements of a matrix C as Cij = 



<yi, coj>; 



calculating the elements of a matrix S as Sij = 



a)j*>; and 



30 



calculating the period matrix R as the solution of 



CR 



= S 



where R satisfies R^ = -I, where I is the identity matrix. 
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32. A method for compressing a mesh representation of a surface, 

the method comprising the steps of: 

conformally mapping the mesh representation to a canonical 

shapes- 
representing the surface position vector on the canonical 

shape as a vector valued function; 

finding the eigen functions . of the Laplacian given by 



decomposing the vector valued function using the eigen 
functions of the Laplacian; 

removing the high frequency components of the decomposed 
vector valued function; and 

storing the low frequency components of the decomposed 
vector valued function as a compressed image of the mesh 
representation, M, of the surface- 

33. The method of claim 31, wherein the surface is a genus-zero 
surface and the canonical shape is a sphere and wherein the step 
of conformally mapping the genus-zero surface to the sphere 
includes the following steps: 

a) computing the Gauss map, mapping M to S^; 

b) computing the Laplacian at each vertex u of M, A(t)(u); 

c) projecting A<t)(u) to the tangent space of <|){u) e S^; 

d) updating (|>(u) along the negative projected A(|)(u); 

e) computing the center of mass of X(()(u), mc(<()); 

f) shifting the center of mass to the center of S^; 

g) renormalizing (|)(u) to be on S^; and 

h) repeating steps b)-g) for all vertices, until the 
projected Laplacian equals zero. 
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34. The method of claim 32, wherein the Laplacian is 



35. A method for compressing a mesh representation, M, of a 
5 surface, the method comprising: 

receiving. the mesh representation; 

conformally mapping mesh representation to a canonical 

parameter domain forming a mapped surface.; 

computing a conformal parameterization of said mapped 
10 surface; 

computing a level set of Gaussian curvature and mean 
curvature for said mapped surface, wherein said level set of 
Gaussian curvature and mean curvature is a function of said 
conformal parameterization; and 
15 storing said conformal parameterization of said mesh 

representation and said level set corresponding to said mesh 
representation. 

36. The method of claim 34 wherein the step of computing said 
20 conformal parameterization includes the step of: 

computing a homology basis {yi, ... , yzq] of said mesh 

representation, wherein computing said homology basis includes 
the steps of: 

computing a set of first and second boundary matrices 
25 di, 82 for said mesh representation; 

forming a matrix D, wherein D = 3292^ + di^di; 

computing the Smith normal for said matrix D; and 

computing the eigenvectors of D corresponding to zero 
eigenvalues that are equal to {yi, 72, ... , 729} wherein {yi, y2, ... , 
30 y2g} is a homology basis of said mesh representation. 
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37. The method of claim 35 wherein the step of. computing said 
conformal parameterization further includes the step of computing 
a harmonic 1-form basis, wherein the step of computing said 
harmonic 1-form basis includes the steps of: 

5 calculating the value of Ci^ = -yi • Yjr i/ j 2, 

... , 2g wherein yi, yj is the homology basis of the mesh 
representation respectively; and 

solving the linear system Sui = 0, Acoi = 0, <ci)i, .yj> = .-yi • yj 
for CDi, wherein {coi, C02f - f ^2q) is the desired harmonic 1-form 
10 basis of the mesh representation. 

38, The method of claim 36 wherein the step of computing said 
conformal parameterization further includes the step of computing 
the fundamental domain Dm of the mesh representation, wherein the 

15 step of computing the fundamental domain Dm of the mesh 
representation includes the steps of: 

selecting an arbitrary face, fo e M, and set Dm to fo; 
setting 30^ = 3fo; 

placing all neighboring faces of fo that share an edge with 
20 fo in a queue, Q; 

while Q is not empty: 

reading the first face in Q; 
setting 5f = eo + ei + e2; 
setting Dm = Dm U f; 
25 finding the first ei e dt such that -ei e SDm; 

replacing -ei in SDMr respectively, by {©(i+i), 

e(i+2) }; 

putting all the neighboring faces that share an edge 
with f in the mesh representation and that are not in Dm or Q into 
30 Q; and 
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removing all adjacent oriented edges in dD^ that are opposite 
each other in sign. 

39. A method for remeshing a mesh representation, M, of a 
5 surface, the method comprising the steps of: 

a) subdividing the mesh, using loop subdivision; 

b) simplifying the mesh using edge collapse that uses a 
minimum edge length criteria. 

c) repeat steps a) and b) until all angles on the mesh 
10 representation, M, are acute; and 

d) outputting the remeshed representation, M, of the 
surface. 

40. A method for converting a mesh representation, M, of a 
15 surface, the method comprising the steps of: 

receiving the mesh representation, M; 

computing a conformal representation of the mesh 

representation, M; 

computing the gradient of the mesh representation, M; 
20 computing the zero points of the mesh representation, M; 

decomposing the mesh representation, M, into canonical 
patches using integration lines along the gradient and through 
the zero points; 

conformally mapping the canonical patches onto a respective 
25 rectangle; 

constructing a tensor product spline surface on the surface 
of the rectangle; and 

matching the control points on the boundary to make the 
parameterization globally smooth . 
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41. The method of claim 39 wherein the mesh representation, M, 
is discontinuous on the boundaries and wherein the 
parameterization is continuous on the boundaries. 

42. A method for conformally mapping a mesh representation of a 
scanned medical image, M, of a body part to a canonical surface, 
the method comprising the steps of: 

a) computing the Gauss map, mapping M to S^; 
* b) computing, the Laplacian at each vertex u of M, A(j){u); 

c) projecting A^{u) to the tangent space of <|)(u) e S^; 

d) updating <t){u) along the negative projected A4»(u); 

e) computing the center of mass of (}){u), mc ((])) ; 

f ) shifting the center of mass to the center of S^; 

g) renormalizing (|>(u) to be on S^; and 

h) repeating steps b)-g) for all vertices, until the 
projected Laplacian equals zero. 

43. A method for animating a mesh representation, .M, of a 
surface from a first mesh representation. Mi, to a second mesh 

20 representation, M2, the method* comprising the steps of: 

removing at least one feature point common to both Mi and M2; 
doubling each of Mi and M2; 

conformally mapping Mi and M2 to first and second canonical 
surfaces, forming first and second mapped surfaces; 
25 computing a holomorphic 1-form for said first and second 

mapped surfaces; 

computing a cohomology basis of said first and second mapped 
surfaces; 

locating first and second zero points on said first and 
30 second mapped surfaces; 
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computing first and second gradients of said first and 
second mapped surfaces at said first and second zero points, 
respectively; 

decomposing the first and second mesh representations^ Mi and 
5 M2/ into first and second sets of canonical patches using 
integration lines along the first and second gradients and 
through the first and second zero points, respectively; 

conformally mapping the first and second sets of patches 
onto first and second rectangles, respectively; 
10 matching said first and second mapped patches on said first 

and second rectangles to form a map therebetween; 

selecting at least one control point on each of said first 
and second mapped patches on said first and second rectangles; 
and 

15 using a BSpline to generate a smooth transition of said 

mapping. 

44. A method of texture mapping on a mesh representation, M, 
representative of a surface, the method comprising the steps of: 
20 removing at least one feature point on M; 

doubling M; 

conformally mapping M to a canonical surface, forming a 
mapped surface; 

computing a holomorphic 1-form for said mapped surface; 
25 computing a cohomology basis of said mapped surface; 

locating zero points on said mapped surface; 
computing a gradient of said mapped surface; 

decomposing said mesh representation, M, into at least two 
canonical patches using integration lines along said gradient and 
30 through said zero point, respectively; 

conformally mapping said at least two canonical patches onto 
at least two rectangles; 
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growing said at least two rectangles until respective . j 
boundaries between said at least two rectangles meet; -J/ 
fixing the boundaries between said at least two recta^^J.es; 

and 

5 solving the Dirichlete problem for the uncovered regions of 

said at least two rectangles. 

45- A method of volumetric harmonic mapping of a mesh 
representation, M, -of a 3-dimensional manifold, the method ^ 
10 comprising the steps of: 

computing the harmonic energy of a mapping f : M-^ R^, where 

the harmonic energy is given by E(f)= X^uv||f(^)""f(v)|| 

[tt^M 

and =-^^QCot(6) , where 9 is the dihedral angle opposite to the 

given edges; 

15 minimizing the harmonic energy using the conjugate gradient 

method to obtain the harmonic mapping, f . 

46. The method of claim 44 wherein the representation M of the 
3-dimensional manifold is a magnetic resonance image of a body 

20 area of interest, and wherein the mapping f : is a map onto 

a canonical sphere of the body area of interest. 

47. The method of claim 45 further including using said map onto 
said canonical sphere of said body area of interest for surgical 

25 planning. 
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